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Measurements of b and c quark asymmetries using data collected at LEP 1 are described. The relative 
merits of each of the individual techniques used is emphasised as is the most profitable way of combining 
^ ' them. Effects of radiative corrections are discussed, together with the impact of these measurements on 

D , global electroweak fits used to estimate the expected mass of the Higgs boson. 

(D 
> 

X 
S3 



1 Motivation for Measuring Quark Asymmetries 



The accurate determination of the forward-backward asymmetries, ^4fb, of quarks serves to test the structure of 
Standard Model (SM) |Q couplings to fermions. They also probe radiative corrections to the SM and consequently 
allow greater precision when predicting unknown parameters of the model. These are increasingly used to 
constrain uncertainties on the mass of the Higgs bosonQ. 

Global fits to electroweak data assume the SM structure of Z couplings to leptons (e, fir) and both up (u, c) 
and down-type (d, s, b) quarks. Given recent, highly accurate, lepton measurements from the r polarisation^) 
and purely leptonic forward-backward asymmetries]^], a similar precision in the quark sector is needed to confirm 
the internal consistency of the model. 

The suite of complimentary measurements described here provide such a precision for both up and down- type 
quark families. Performing these measurements at LEP 1 offers several advantages. The sensitivity of initial 
state couplings to the effective weak mixing angle, sin 2 6^ , is compounded by large, measurable asymmetries 
from quark final states close to the Z. Heavy quark asymmetries in particular are especially favourable, as flavour 
and direction of the final state quark can be tagged with greater ease than is the case with lighter quarks. 

2 Definitions and Experimental Issues 

In the SM, the differential cross-section for the process e + e~ — > / + /~ can be written as : 

where A^ B defined to be the forward-backward asymmetry for fermion flavour, /. It can be expressed as : 

4 B =- A A e A f (2) 



where Aj is the polarisation of the fermion concerned : 



A, = -^— 2 =1-H (sin 2 ^ + C f ) (3) 
1 + x JJ 



where x is the ratio of the vector and axial couplings of the fermion to the Z. This final form, separates the 
terms containing sensitivity to parameters of the SM, such as (rnt, toh) through sin 2 #^, from vertex corrections 
in Cf. The latter are typically of the order of ~ 1% for b quarks B. 

For hadronic decays of quarks, the precise direction of the final state fermion is not accessible experimentally, 
and so the direction of the thrust axis is usually signed according to methods of correlating the charge of the 
quark its final decay products. Asymmetries are of the order of ~ 10% for b and c quarks but are diluted by 
several effects. These are caused primarily by the correlation method mistagging the quark charge, or by B°B° 
mixing or cancellations between other quark backgrounds. 

Consequently, the methods of measurement described here represent different compromises between rates of 
charge mistag, and the efficiency+purity of the flavour tagging procedure. With the increasing sophistication of 
analyses, several methods of tagging the charge and flavour of decaying quarks are available. These are applied 
either singly or in combinations. 

Minor complications arise when interpretating the results of these analyses, in the form of corrections to pure 
electroweak predictions of heavy quark asymmetries. For example, quark mass corrections to the electroweak 
process are generally small, eg. representing shifts of 0.05% in the calculation of Ap B , and are well understood 
theoretically Larger and more problematic corrections arise from hard gluon emission in the final state. 
These so-called QCD corrections are mass, flavour and analysis dependent and so their treatment and associated 
uncertainties must be handled with care. 



3 Latest Techniques for Measuring Ap B 
3.1 Semileptonic Decays of Heavy Quarks 

The "classical" method of measuring heavy quark asymmetries relies on differences in the momentum and 
transverse momentum, (p, pt), spectra of leptons arising from semileptonic b and c quark decays. This method 
benefits from an unambiguous charge and flavour tag in the case of unmixed b and c hadrons. It however suffers 
from several disadvantages when used on a more general sample. 



The methods reliance on a pair of correlated inputs, such as lepton (p, px) a lone, leads to a dependence of 
the method on the precise values of the branching ratios, BR(6 — > I) and BR(fe — > c — > I), semileptonic decay 
modelling and c branching fractions assumed in Monte Carlo simulations. Some uncertainties are constrained by 
full use of measurements from lower energy experiments, according to the prescription detailed in[Q|. However, 
the extrapolation of such measurements to LEP energies, and into the different environment of a fragmented 
jet containing a 6-hadron, leaves residual systematic uncertainties. As the purity of the sample, and its charge 
mistag, are determined by Monte Carlo simulation alone these residual uncertainties are correlated between the 
4 LEP experiments. 

The relatively small value of the (b — > I) branching fraction means that an overall b tagging efficiency of the 
order of ~ 13% is typically achieved with reasonable purities of eg. ~ 80% in the case of||. Effects of B°B° 
mixing, "cascade" b decays (b — > c — > I), backgrounds from c decays (c — > I) and other light quark sources, 
including detector misidentification, all serve to dilute the otherwise excellent mistag rate. 

The systematic impact of such effects can be severely reduced by fitting simultaneously for the the time- integrated 
mixing parameter, jj-PL together with the b and c asymmetries in a "global" fit. Otherwise both \ an d A FB are 
inputs to the determination of Ap B . Their values and uncertainties are either fixed by experiment or, in the case 
of A FB , set to their SM expectation and the dependence on A FB used as input to subsequent electroweak (EW) 
fits Q . Uncertainties can be further minimised by adding input information which discriminates between b and 
c events. In such analyses || Q information coming from lifetime tags, using silicon vertex detectors (VDET's), 
and event shapes amplify the discrimination from semileptonic (p, px) spectra. 

The LEP experiments make use of such additional inputs, and extra fit quantities, to varying degrees as is 
summarised in Table |l|. Analyses which make use of the classical semileptonic (p, px)spectra, and complement 



Inputs 


lepton [p, px) 


lifetime 


jet shapes \ 


Expt. 


ADLO 


DO 


DO DL 


Outputs 




X 




Expt. 


ADLO 


AO 


DO 



Table 1: Additional inputs and fit parameters for the semileptonic decay measurements of Ap B from the 4 LEP 
collaborations ALEPH, DELPHI, L3 and OPAL (ADLO). 



it with lifetime tags in this way, are the most powerful, statistically and with greater systematic control. For 
example, the OPAL measurement^ uses lepton (p, px)and event-shape information as inputs to a neural net 
(NN) b tagging algorithm. In addition, it utilises a largely orthogonal c tag, based on lifetime information of 
jets in the event, combined with the impact parameter significance and detector identification criteria of the 
lepton. The output of these two neural nets is shown in Figure [j] where the strong separation between sources 
of leptons in hadronic events is clearly evident. The separation between b and c lepton sources, and the more 
limited distinction between those and other background sources, enables both an precise determination of Ap B , 
X and the most accurate measurement of ^4 FB from the same sample of events. The net gain of such a method 
is an approximate ~ 25% improvement in the statistical sensitivity of A FB and a ~ 25% improvement in that of 
A FB . A summary of the current results for A FB from semileptonic measurements at LEP is given in Table || 



Semileptonic Measurements of A FB 


Experiment 


A FB stat. syst. total 


ALEPH 

DELPHI 

L3 

OPAL 


0.0965 ±0.0044 ±0.0026 ±0.0051 
0.0979 ±0.0065 ±0.0029 ±0.0071 
0.0963 ±0.0065 ±0.0035 ±0.0074 
0.0910 ±0.0044 ±0.0020 ±0.0048 


Lifetime and Jetcharge Measurements of A FB 


Experiment 


Ap B stat. syst. total 


ALEPH 

DELPHI 

L3 

OPAL 


0.1040 ±0.0040 ±0.0032 ±0.0051 
0.0979 ±0.0047 ±0.0021 ±0.0051 
0.0855 ±0.0118 ±0.0056 ±0.0131 
0.1004 ±0.0052 ±0.0044 ±0.0068 



Table 2: Summary of latest measurements of A FB from the 4 LEP experiments. 



^ X is defined to be the probability that a B meson has oscillated to a B meson by the time of its decay. 
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Figure 1: Two-dimensional summary plot of the c tag neural net output vs. that of the b tag NN output for each 
of the expected lepton sources in the OPAL semileptonic asymmetry analysis. The areas of displayed pixels is 
proportional to the bin contents. 



3.2 Lifetime Tagging and Jetcharge Measurements 

An alternative, complementary technique to measure Ap B is based on lifetime information from silicon vertex 
detectors. This is then combined with a fully inclusive charge correlation method, referred to as the "jetcharge" 
technique. This method was initially pioneered using samples of untagged hadronic events containing all types 
of quark flavours accessible at these energies Q. As a consequence of the low semileptonic branching ratios, 
such inclusive measurements are almost entirely uncorrelated from semileptonic measurements and so can either 
be combined or used as a cross-check for consistency of measurements between different methods. 

The jetcharge method is based upon the correlation between leading particles in a jet, with that of the 
parent quark. A hemisphere based jetcharge estimator is formed using a summation over particle charges, q, 
weighted by their momentum, p : 

t . ^"' f f' , (4, 
E f-™|fl.f |« 

and analogously for Qb- The k parameter is used to optimise the measurement sensitivity. The charge flow 
between hemispheres, namely Qfb = Qf — Qb is then used to sign the direction of the thrust axis. Currently 
all LEP collaborations use this method, and the above formalism. 

The method benefits from many of the systematic studies performed in untagged samples [^8| , especially to 
understand the degree of charge correlation between hemispheres in background events and their light quark 
parents. The recent Ap B analysis carried out by DELPHI also illustrates this method [fL0|. In addition to 
jetcharge information, OPAL makes uses of a weighted vertex charge method[^2|. This quantity is a weighted 
sum of the charges of tracks in a jet which contains a tagged secondary vertex, ie : 

<2Vtx = ^ t qi (5) 

tracks— i 

where g, is the charge of each track, i and is related to the probability that the tracks comes from the 
secondary vertex relative to that it came from the primary. The latter probabilities are determined using 
impact parameter, momentum and multiplicity information. An estimate of the accuracy of the q vtx charge 
estimator is derived from its variance. 

Selecting hemispheres with | q vtK \> 1.4 x a q + 0.2 largely removes neutral B° mesons, and those with poorly 
measured vertex charges. This also leads to a severe reduction in the size of the event sample, leaving only 
~ 13, 000 events out of a total, untagged input of roughly 4 million hadronic events. Hence, the contribution of 
the vertex charge measurement, when combined with the jetcharge determination of ^4p B , is relatively low. 



More significant improvements in statistical precision and control of systematic uncertainties can be obtained 
from the variety of new techniques summarised in Table ^. Experiments make use of these techniques to varying 



Expt. 


Improvement 


A L 


Fit data as a function of tag purity. 


A L 


Increase acceptance beyond central region. 


ADLO 


Fit asymmetry as a function of angle. 


A 


Fit data using range of k values. 


ADO 


Extract b quark mistag factor from data. 


AD 


Extract lighter quark mistags from data. 



Table 3: Recent improvements to the lifetime and jetcharge method for measuring Ap B from the 4 LEP collabo- 
rations ALEPH, DELPHI, 13 and OPAL (ADLO). 



degrees, the most significant of which being the improvement in statistical sensitivity gained by fitting to the 
asymmetry as a function of angle. With increased statistical precision, comes the need for improved systematic 
control. The most important of these being the extraction of the charge mistag factor for b quark from data. 
ALEPH, DELPHI and OPAL now perform this extraction while ALEPH also extracts it as a function of the polar 
angle of the thrust axis. This takes into account particle losses close to the edge of detector acceptance. 

The output distributions from the ALEPH measurement are shown in Figure pi The asymmetry measurement 
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Figure 2: Extracted fit variables from the ALEPH lifetime+jetcharge measurement of Ap B . Columns represent 
the b charge mistag factor and sample b purity as a function of angle in samples with increasing b purity. Filled 
and open circles represent n values of 0.5 and 2.0 in equation (4) respectively. 

is made separately in each bin of polar angle, n and b sample purity before being combined. Some points of 
interest in such new analyses include the increased statistical power arising from the bins at low angles, even 
those where the thrust axis lies outside the VDET acceptance. The large value of the asymmetry in these regions 
compensating for the low tagging efficiencies. The increase in sample b purity at large cos 8 for high b purity 
samples, is due to the loss of tracks at the edge of the VDET acceptance. This affects b events the least, as 
more tracks with large pt to the thrust axis continue to tag the event. This is true to a much lesser extent for 
lighter quark flavours. 



Measurements of the b asymmetry using the lifetime and jetcharge method are also summarised in Table 0. 
It is interesting to note that such measurements, whilst providing similar sensitivity to A FB as semileptonic b 
decays, as yet do not provide the possibility of analogous measurements of Ap B . 



4 Latest Techniques for Measuring the c Asymmetry 

In contrast to the incremental progress in the field of b asymmetries, measurements of the corresponding quantity 
in c decays have improved dramatically in recent years. Of the 4 LEP collaborations, DELPHI, L3 and OPAL 
have determined the c asymmetry as an output of global semileptonic fits to b and c decays^, |^, |?J whereas 
ALEPH, DELPHI and OPAL have performed the same measurement using fully reconstructed samples of D meson 
decays [|l3[ As semileptonic measurements are discussed in Section 3.1, only the latter are described 

here. 

The method of exclusively reconstructing D decays aim to use as many channels as possible by reconstruct- 
ing the D° through its decay to K~tt + . The D° is generally reconstructed Jl3[ by taking all 2 and 4 track 
combinations and a ir° candidate with zero total charge. Those combinations with odd charges are then used 
to form possible D + candidates. Each experiment reconstructs a different subset out of a total of 9 different 
decay modes. The dominant channels however are the D + — > AT~7r + 7r + and D* + modes, with all modes offering 
some statistical power. Similarly, each experiment has different selection criteria in each mode, depending on 
the momentum and particle identification resolutions of the detectors. 

These differences lead to widely varying efficiencies and signal-to-background ratios. For example, the 
DELPHI mass difference distributions for 4 of the 8 selected modes are shown in Figure |3[ An important 
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Figure 3: Difference distributions between the mass of the D* + and the D° candidate for different decay modes. 

advantage of such measurements is the ability to determine the background asymmetry, Apg S , from data, 
mode-by-mode using information from sidebands. 

A major difficulty is encountered when trying to correct for the substantial fraction of events due to B° — > D 
decays for B° — B° mixing. Each D mode is corrected using an "effective" x depending on the expected fractions 
of and B® decays contributing to the mode concerned. These effective factors are determined fromjl6| using 
Monte Carlo simulation and so give rise to systematic uncertainties. The observed asymmetry, A^, is then 



found using : 

^FB = /sig/c^FB + /slg(l — /c)^FB + (1 — /sig)^4 FB S (6) 

where / s i g is the fraction of signal of signal+background events and f c is the fraction of events containing a true 
D meson which are due to c quark events. As far as possible, the sample c purities, f c , are determined from 
data using lifetime, mass and event-shape information in both hemispheres of events containing a D tag. 

Each experiment makes use of lifetime information, with both DELPHI and OPAL using the D momentum 
and jet-shape information respectively in addition, so as to disentangle the substantial contamination from b 
decays. The DELPHI experiment does so in the context of a simultaneous fit to both b and c asymmetries. 
However, besides constraining systematic uncertainties, the precision available of ^4p B is negligible compared to 
methods discussed previously. 

The c asymmetry measurements described here are summarised in Table ^ and indicate that, despite com- 
plex systematics, such measurements remain primarily limited by the low efficiencies and purities of the D 
meson reconstruction. Systematic errors vary widely, with the time-dependence of the background asymmetry 



Semileptonic Measurements of A^ B 


Experiment 


Ap B stat. syst. total 


DELPHI 
L3 

OPAL 


0.0770 ±0.0113 ±0.0071 ±0.0133 
0.0784 ±0.0370 ±0.0250 ±0.0446 
0.0595 ±0.0059 ±0.0053 ±0.0079 


Exclusive D Tag Measurements of A^ B 


Experiment 


Ap B stat. syst. total 


ALEPH 

DELPHI 

OPAL 


0.0630 ±0.0090 ±0.0030 ±0.0095 
0.0658 ±0.0093 ±0.0042 ±0.0102 
0.0630 ±0.0120 ±0.0055 ±0.0132 



Table 4: Summary of latest measurements of Ap B from the 4 LEP experiments. 



remaining merely one of many dominant sources depending on decay mode and experiment. 

5 Radiative Corrections to Asymmetry Measurements 

Several small corrections must be made to the b and c asymmetries extracted from the described analyses. In 
the case of Ap B , QED corrections for ISR and FSR are relatively minor, amounting to -0.0041 and -0.00002 
respectively. Similarly, corrections for pure 7 exchange and 7 — Z interference diagrams give rise to a correction 
of +0.0003. Such corrections, are general in nature and so apply equally to all analyses. 

A more difficult set of corrections involves those needed to correct for the presence of hard gluon radiation 
which can distort the angular distribution of the final state quarks when compared with the pure electroweak 
process. Estimates of such corrections to heavy quark asymmetries have been computed to first and second 
order in a s both numerically |P^] and, most recently, analytically p9[ in different scenarios for either c, b or 
massless quarks. 

A common procedure for correcting and ascribing systematic uncertainties for the LEP heavy quark asym- 
metries has been developed ^o). The more recent analytical calculations indicate several discrepancies when 
compared with the numerical results. These remain to be resolved. Current systematic uncertainties are de- 
termined using a procedure of comparing the effects between first and second order in QCD and by switching 
between massless quarks, and assumptions for the c and b quark masses. 

Further difficulties arise when considering the application of such corrections to individual analyses. Theo- 
retical calculations are typically based on the direction of the outgoing quark, whereas the analyses described 
here use the thrust direction. Further, the sensitivity to hard gluon radiation of data, containing either a 
lepton of a given (p, pt), a reconstructed D meson or purely inclusive events, varies dramatically. Effects of 
non-perturbative QCD and higher-order effects during hadronisation must also be evaluated. The latter ren- 
der QCD corrections both detector and analyses dependent, eg. event shape selections implying an implicit 
dependence on the strength of gluon emission. The correction to be applied to a given analysis is derived from : 

^FB = (1 - Cb,cSb.c) I no QCD (7) 

where C{, jC represents the QCD correction at parton-through-to-hadron level, and Sb, c is the analysis dependent 
modification. Examples of the magnitude of the QCD corrections at the theoretical and experimental levels 
are shown in Table o for the cases of Ap B , determined using semileptonic and lifetime+jetcharge analyses p0[ . 





Lepton Analyses 


Lifetime+Jetcharge 


St min. 
max. 


0.52 ±0.06 
0.74 ±0.07 


0.24 ±0.46 
0.36 ±0.32 


Cb min. 
max. 


1.54 ±0.28 
2.19 ±0.37 


0.71 ±1.36 
1.07 ±0.96 



Table 5: Summary of QCD corrections to Ap B for the different analysis methods. 



The constants are evaluated in terms of Monte Carlo simulations, before and after experimental cuts. The 
hadronisation dependence of corrections is included as a systematic uncertainty by comparing results from both 
the HERWIG[^2] and JETSET|2^] models. It is seen from comparing parton and hadron level that the effect of 
hadronisation is to reduce the magnitude of the QCD correction]^]. 

It is important to note that the corrections for lifetime±jetcharge measurements are negligible. Significant 
corrections are observed for both semilcptonic,and D tag measurements of both A FB and Ap B . Jetcharge mea- 
surements are immune to such corrections as the 6-quark charge mistag factor is defined using Monte Carlo with 
respect to the original 66 quark pair orientation, prior to gluon or final state photon radiation, parton shower, 
hadronisation and B°B° mixing. All these effects are therefore included, by construction, in the analyses, as 
far as they are properly modelled in the JETSEt|2^] hadronisation model. 



6 Conclusion and Perspectives 

With the completion of LEP data-taking at energies close to the Z resonance in 1995, the 4 experiments 
(ALEPH, DELPHI, L3 and OPAL), have accumulated large samples of hadronic events. From this data, the 
forward-backward asymmetry of the 6 quark has emerged as the most sensitive single test of the SM at LEP. 
The complementary c asymmetry measurements offer additional precison and a new window on couplings in 
the down-type quark family. The precision from semileptonic measurements of j4 FB is now matched by that of 
lifetime±jetcharge measurements. The electroweak sensitivity of Admeasurements now equals that obtained 
from combined quark asymmetries measured in untagged samples, highlighting the benficial effect of flavour 
tagging. 

However, in light of these measurements great sensitivity to the couplings of the SM, the continuing dis- 
crepancy between electroweak results from LEP and SLD|2^] make it essential to understand whether it is due 
to statistical fluctuations or systematic effects. Separating LEP measurements of ^4 FB into those from the two 
dominant techniques, and conservatively ignoring correlated systematic uncertainties, indicates that there is 
at most a 1.2ct discrepancy between semileptonic and lifetime+jetcharge measurements. This is insufficient to 
explain the LEP-SLD discrepancy but indicates that care must be taken when considering common systematics 
in leptonic decay modelling and fragmentation uncertainties. Further improvements in both 6 and c asymme- 
tries are possible, as both sets of measurements are still dominated by statistics. For semileptonic analyses, the 
benefits of using both lifetime and lepton information are emphasised. In the case of the lifetime and jetcharge 
method, these are most likely to come in the form of improved 6 tagging efficiencies and extensions of tagging to 
lower angles. The situation for improvments to measurements of A FB is more difficult as the number of available 
modes is exhausted, and efficient methods of tagging c events remain to be discovered. 

At this point, without the prospect of significant, further LEP data-taking at the Z, it is important to 
focus upon the latest techniques which offer the greatest sensitivity to the couplings of the SM combined with 
systematic control. The measurements described here obtain combined precisions on the 6 and c asymmetries 
of 2.2% and 7.1% respectively. Hence, the goal of acheiving similar precision on the Z couplings to quarks, as 
that obtained for leptons, has been reached. 
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